Gauge theory of topological phases of matter 
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We study the response of quantum many-body systems to coupling some of their degrees of freedom to small, 
slowly varying external gauge fields, which arise when global symmetries are gauged. This analysis leads 
to a "gauge theory of states of matter" generalizing the well known Landau theory of order parameters. We 
illustrate the power of our approach by deriving and interpreting the gauge-invariant (local) effective actions of 
superconductors, 2D electron gases exhibiting the quantized Hall- and spin-Hall effect, 3D topological insulators 
and axion electrodynamics. 
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During the past few years, various novel states of con- 
densed matter protected by topological properties (but not de- 
scribable by local order parameters) have been predicted and 
found in experiments; see Refs. |l-9] and references given 
therein. The existence of several of these states was predicted, 
conjecturally, in the 90's, based on what might be called a 
"gauge theory of states of matter" 14], but these results 
do not appear to have been widely noticed in the condensed 
matter community; (for earlier work, see also Refs. lll5l42lll ). 
The purpose of this letter is to recall some key elements of the 
theory and to show how it can be used to recover various re- 
cent theoretical predictions and speculations and to interpret 
experimental findings. We limit our analysis to electron gases, 
but the theory also applies to other systems such as cold atom 
gases [11], the primordial plasma in the universe ll22ll . etc. 

The key idea is to explore states of electron gases by an- 
alyzing their response to turning on external gauge fields 
and determining their effective action or effective free en- 
ergy. The Pauli equation of a non-relativistic, spinning elec- 
tron reads ihdt'^t — Ht'^t, where the Hamiltonian is given 
by Ht = 2^11^]^ + ^(2^7 ^ is a potential and *t(a;) a 
two-component spinor. This equation is invariant under global 
phase- and SU (2)-transformations of the state function ^"4. In 
accordance with Noether's theorem, there are two conserved 
current densities, the electric- and the spin current density. We 
promote these global symmetries to local gauge symmetries 
by introducing U{1)- and S'C/(2)-gauge fields (vector poten- 
tials), a and w. These gauge fields describe, among other 
things, effects of external electromagnetic fields, spin-orbit 
and Zeeman interactions, effects due to the curvature of the 
sample (as well as dislocations and disclinations), and the 
influence on electronic properties of the motion of the ionic 
background harboring the electron gas. In order to take into 
account the latter, the Pauli equation is written in moving co- 
ordinates corresponding to a (divergence-free) velocity field 
V, which yields contributions to the gauge fields a and w, as 
described below. The resulting Pauli Hamiltonian is given by 
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+ ao + WQaK- (1) 



Here, ctx, K ~ 1, 2, 3, are the usual Pauli matrices; a{x, t) — 
{x, t) + f t) + mv{x, t), where A'^°^ may describe 



a static, homogeneous external electromagnetic field and geo- 
metrical properties (curvature) of the sample, A'^™ is the vec- 
tor potential of fluctuations in the electromagnetic field, and v 
is the velocity field describing the motion of the ionic back- 
ground of the electron gas. (That the velocity field appears as 
a contribution to the U (l)-vector potential can be understood 
by recalling that the Lorentz- and the Coriolis force have the 
same general form.) The constant e is the elementary elec- 
tric charge, c the speed of light, and m the electron mass. 
The time-component is ao(a;,t) — ^Lp{x,t) — ^'iP{x,t) + 
j;p{x, t), where ip is the electrostatic potential, p the pressure 
and p the density of the ionic background. Furthermore, the 
5'C/(2)-gauge field, w, is given by 
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w^aK^jr{B + —^Av)-'^ + Wo-<J, (3) 
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with n = - 2^)(|-^ + ^^)- The first term in EqJZ) 

describes spin-orbit interactions and Thomas precession 12311 . 
The fields Wq and describe magnetic exchange interac- 
tions (Wq is the "Weiss exchange field"), while fl is the spin 
connection describing parallel transport of spins in a curved 
background and interactions with dislocations and discUna- 
tions. The fields E and B are the electric field and the mag- 
netic induction, respectively. We observe that v is gauge- 
invariant under C/(l)-gauge transformations and that H and 
W transform homogeneously under SU (2)-gauge transforma- 
tions, while (l transforms inhomogeneously. All vector quan- 
tities in Eqs. (2) and (3) are expressed in the local coordinate 
systems. 



Introducing the covariant derivatives Dj 

iao + iwQan 

Pauli equation in the form 
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iw^ax and Dq = hdt + iao + iw^aK^ we can write the 
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where g'^ is the metric tensor of the sample background and 
g its determinant (see Ref. ifioll ). This equation displays full 
U{l)sm X SU{2)spin gauge invariance. It turns out to be the 
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Euler equation corresponding to the action 

Soi¥,^;a,w)= [ dt [ ^dxhlix) ■ Do-i'tix) 



(5) 



In order to describe electron-electron interactions, we add 



dt / ^dx\^t{x)\^U{x-y)\^tiy)\^ 



where W is a two-body potential. The total action functional 
is given by S" = + ^int- 

In what follows, we propose to study the response of an 
electron gas with a bulk mobility gap (an insulator) to turning 
on external gauge fields a and w. (More generally, one may 
also gauge "emergent" symmetries of such systems.) For sim- 
plicity, we only study ground state properties of such systems 
in this letter Thus, we determine the form of the effective ac- 
tion at temperatures T « 0. This can be done by considering 
the expectation of the propagator, Ua,w{t, s), of an interact- 
ing electron gas (with two-body potential U) from time s to 
time t, in the presence of time-dependent gauge fields, in the 
ground state, |i^o), of the gas. One defines a "partition func- 
tion" Z{a,w) := \iT[is^-oo{ipa\Ua,w{t, s)\lpq). This quantity 

can also be expressed as a functional integral by promoting 
{x) and ^'J (x) to Grassmann variables and performing the 

Berezin integral Z{a, w) = const • / D'i'' D'^/e^:^^'^'* •'^•"•'"l 
The effective action is defined as 



Seff{a,w) = —ihlnZ{a,w). 



(6) 



It has the following general properties: 

(1) It is the generating function of Green functions of the 
electric current density and the spin current density s^: 



dSeffia,w) 

da^{x) 
9S.n{a, w) _ 



(7) 
(8) 



while higher derivatives yield connected Green functions of 
these current densities. (The derivative of S'eff with respect to 
the metric gij is the expectation value of the stress tensor) 
(2 j It is gauge invariant: 



where x is a real-valued function and U denotes a space-time 
dependent rotation in spin space. We note that electromag- 
netic gauge invariance and electric current conservation are 
equivalent, while iS'J7(2)-gauge invariance is equivalent to the 
property that the spin current is covariantly conserved 1 10]. 

(3) Assuming that connected current Green functions have 
appropriate cluster properties (which is the case for gases with 



a mobility gap above the ground state energy, such as insula- 
tors) and passing to the limit of large distance- and low fre- 
quency scales (i.e., the scaling limit), S'eff (a, w) can be written 
as a sum of integrals over local polynomials in a and w and 
derivatives thereof. These polynomials are gauge-invariant up 
to total derivatives, which are cancelled by appropriately cho- 
sen boundary terms; (this is commonly called "anomaly can- 
cellation"). In a system confined to a region A of space-time 
with non-empty boundary i9A, the effective action thus takes 
the form 

n 

where n denotes the scaling dimension of these terms; (the 
dimensions of a derivative and of a gauge potential are = 1). 

(4) In order to explore physics in the (long-distance-, low- 
frequency) scaling limit, it suffices to retain the leading (most 
relevant) terms in the expansion of S'eff. 

While the starting point of Landau theory is the identifi- 
cation of global symmetries of a system and of a space of 
order parameters on which these symmetries act, along with 
the allowed patterns of spontaneous symmetry breaking (see 
1 24, 2^), the starting point of our "gauge theory of states of 
matter" is the idea that all global (fundamental or emergent) 
symmetries shall be gauged and the response of the system 
to turning on the corresponding external gauge fields be an- 
alyzed. This enables one to identify states of matter not de- 
scribable by local order parameters, such as states exhibiting 
"topological order", and to analyze various transport coeffi- 
cients. In this letter, we illustrate the general theory on the 
special example of electron gases with a bulk mobility gap. 

(2 + 1)D examples, (i) We start by considering a two- 
dimensional electron gas subject to a strong, uniform mag- 
netic field 5^°) = VAA'^™'-'^^ perpendicular to the sample sur- 
face (and defining the local z-axis). The vector potentials w 
and V are set to zero. We assume that, for appropriate choices 
of the external field B^°\ the bulk Hamiltonian of the gas has 
a mobility gap above the ground state energy, with the spins 
of the electrons aligned in the direction of . We propose 
to study the response of the electron gas to small fluctuations 
in the electromagnetic field and in the curvature of the sample 
surface. The total electromagnetic vector potential is denoted 
by A""' = + A. If the sample surface of the gas has 

non-vanishing Gauss curvature, the C/(l)-gauge field con- 
tains a contribution describing parallel transport on the sample 
surface (rotations around the local z-axes) besides A. Equa- 
tion ^ then leads to an effective action of the form 



Seff(a) 
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dtdxe^'^Pa^.d^ap + T{a\dA) 
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+ r(a|aA), (9) 



where an is the Hall conductivity, e^'^'^ the Levi-Civita anti- 
symmetric tensor, Aq ^ Lp the scalar potential, and K is the 
Gauss curvature of the sample. This action describes the 
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Quantum Hall Effect (QHE): 



(10) 



Curvature effects in the context of the QHE have first been 
described in Refs. [Til] (see also Ref. ll26ll ). and have re- 
cently attracted renewed interest (see, e.g., Ref. The 
edge action r(a|aA) is the well known anomalous chiral ac- 
tion in (1 + 1) dimensions, which cancels the gauge anomaly 
of the first (bulk) term (the (2 + l)-dimensional Chern-Simons 
action) in Eq. It is the generating function of Green func- 
tions of chiral electric edge currents propagating along the 
boundary 9A of the sample, see Ref. |28] and [ isi [l9ll . Ef- 
fects of the motion of the sample harboring the electron gas 
can be accounted for by adding the velocity field v describing 
the motion of the sample to the J7(l)-connection ^1^,^11']. 

(ii) Introducing an S'C/(2)-gauge field, w, describing 
exchange- and/or spin-orbit interactions, we find the effective 
action appropriate to describe the Hall effect for the spin cur- 
rent, as described in detail in Refs. |Tdi Ul- In the scaling 
limit, the most general effective action for two-dimensional 
systems is given by 



^eff(fl, W) = (2A2 



dtdx{d 



•T\2 



X J dtdxTr{U + Wf + ^ y dtdxe^'^'Pa^.d^ap 

k f ( 2 

— / dtdxe'^^^Tr iwp_d^Wp + -w^w^Wf 

dtdx{l/^i)B'^ + J dtdxeE"^ 



edge action depending on a|aA and w\gf^. 



(11) 



Here, and in the following, the gauge coupling constants are 
absorbed into the definitions of the magnetic permeability ji 
and the dielectric constant e, aF is the transverse part of a and 
A the London constant of a superconductor; the field wq is 
given by Eq. (|3]l. The coefficients x x proportional 
to susceptibilities. The fourth and fifth term on the right hand 
side are Chern-Simons terms that have gauge anomalies at the 
boundary of the sample and require the addition of an "edge 
action" canceling these anomalies. Of course, for most sys- 
tems, only some of the terms on the right side of Eq. (fTTT i are 
present; (e. g., the first and the fourth term tend to exclude 
one another, because superconductivity requires time reversal 
invariance, while the QHE requires its breaking). 

The coefficient k appearing in front of the second Chern- 
Simons term in Eq. (fTTT i must be an integer [29]. The bound- 
ary action canceling the anomaly of this term is the generating 
function of Green functions of chiral current operators gener- 
ating an SU (2) current algebra at level k. The unitary repre- 
sentations of an SU (2) current algebra at level k are labeled 
by a spin quantum number s=:0,|,...-|,((fc + l) inequiva- 
lent irreducible representations). Level fc = corresponds to 
topological insulators without chiral boundary currents, while 
fc = 1 corresponds to chiral edge spin currents carried by 



quasi-particles of spin i . The range fc > 2 would be relevant 
only if there existed stable quasi-particles at the edge with spin 
■| > 1. Recalling Eq. (|2|, one observes that spin-orbit inter- 
actions can excite chiral edge spin currents. For alternative 
results on 2D topological insulators see also Ref. 1,30.1 . and 
referrences given there. 

Further insights gained from the analysis of Eq. ( fTTT i in- 
clude a general version of the Goldstone theorem for 2D sys- 
tems, a duality between insulators and superconductors in two 
space dimensions, and the prediction of a Hall effect in rapidly 
rotating 2D atom gases (with v playing the role of the C/(l)- 
vector potential); see Ref. 

(3 + \)D examples, (i) Superconductors: We first set the 
SU (2)-gauge field w to 0. The effective action becomes 

^eff(a) = ^ / dtdx{d^f + 32^ / dtdxe'^^'^'F^^Fp,, 

(12) 

where — dpa^ — d^a^ is the electromagnetic field ten- 
sor. The second term in Eq. (fT2] i is a topological term; (it is 
really a surface term). For general values of 7, it breaks par- 
ity and time reversal symmetry. We note, however, that, for 
an infinitely extended sample, and with continuity conditions 
imposed on a at infinity, the range of values of this term is 
given by 7 X Z. Thus, for 7 0, vr, exp[iS'eff(a)] preserves 
parity and time reversal in the bulk. 

Equations (fT2l) and (|7]i lead to the London equation: 
2 

em\T 
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(13) 



where = —mc/e^n, with n the condensate density. 

(ii) Next, we consider insulators, which are materials with 
a bulk gap for which the condensate density n vanishes, and 
determine their effective actions in the presence of C/(l)- and 
S'C/(2)-gauge fields. For simplicity, we set v — 0. The effec- 
tive action then reads 



-Seff (a, w) ^ - 
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less relevant terms. (14) 



Here, {Fyj)py = d^Wy - 9^w^ + ^[wp,Wy\. We assume 
that the sample A has the geometry of a "slab" and denote 
its boundary by dh. From Stokes' theorem we find 

1 



dtdxE B ^ - I dtdxe>"'Papdyap =: 27r'^raA(a), 
A 2 



(15) 



i / dtdxe''''''^Ti[{F^)^,{F^)p,] 

4 JA 



dtdxe'^'^PTi 



OA 



WpdyWp + -w^w„Wp =: 47r^raA(?i)). 

(16) 
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On the right hand side of Eq 

non-abelian Chern-Simons term. 

1 



. ( fTSI l we recognize the 
We note that, because 
32^2 / dtdxTr[£^'''P''{F^)^,^{F^)p„] is an integer, parity in- 
variance of the bulk impHes that = or tt, and conversely. 

The boundary term given by 'nTQ/^{a) is the effective ac- 
tion of a charged, 2-component (massless) Dirac fennion 
llsill . (Quasi-particles with the same properties also appear 
in graphene, see Ref. [15]). Thus, -f/ir determines the number 
of species of charged Dirac fermions propagating along the 
boundary. The term ttTqa (w) is the effective action of a rel- 
ativistic fermion with "SU{2) isospin". We conclude that the 
gauge-invariant action (fT4l i predicts topological insulators, for 
7 and/or 9 0; (j,9 = corresponds to ordinary insulators). 

Among further applications of the general ideas described 
here we mention analyses of the Einstein-deHaas-Barnett ef- 
fect, of vortices in superfluids (see Ref. Ill ill ), or of the phe- 
nomenon of sonoluminescence. 

Axion electrodynamics. If we promote the coupling 7 
in Eq. (fT4l ) to a dynamical variable (field) (f), i.e., replace 
/ dtdxE • / dtdxe^"'P''F^,^Fp„ by 
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dtdxe''''P''{^ + l(t>)F^,Fp,, 



(17) 



where / is a parameter with the dimension of a length, and add 
the term ^ J dtdx ^di^i_4)dP(j) + U{(t)) , then we obtain an ac- 
tion containing a coupling of electrons to an axion field; see 
Refs. lfl3i[l4l] . The axion potential U{(j)) is usually periodic 
in (p with period 2-k/I. If [/ 7^ then one predicts the ex- 
istence of domain walls across which the value of the axion 
field changes by an integer multiple of 27r/Z . These domain 
walls, which, for entropic reasons, must occur in the bulk of 
axionic topological insulators (or axionic topological super- 
conductors), support massless charged modes. 

(4 + 1)D examples. The QHE has a (4 + 1) -dimensional 
cousin first studied in Refs. ifo, [l4ll . Let us consider a five- 
dimensional system confined to a slab Q < x'^ < L consist- 
ing of very heavy charged four-component Dirac fermions. If 
these Dirac fermions are coupled to an external electromag- 
netic vector potential A*^^) and are then integrated out, the ef- 
fective action, as given by Eq. becomes 



with S 



fMiA^'^) = '^!dtdx{F(^)Y^{F^^)),, the (4 + 
1) -dimensional analogue of the Maxwell action (F^^'' = 
d^,A^^^ - d^A'^^\ and S'^^|(A(5)) proportional to the five- 
dimensional Chern-Simons action 



N 
967r2 



dtdxeP^^P'A^^^FlfFll\ 



(18) 



where N ~ 1,2,... is the number of fermion species. The 
boundary term rar((^^^-* |ar) must be introduced in order 
to ensure the gauge invariance of the effective action in the 
slab geometry. It describes massless chiral fermions on the 



(3 + 1) -dimensional "top" and "bottom" boundary compo- 
nents (or "branes") of the slab. These chiral fermions may 
acquire a mass through tunneling between the two boundary 
components. 

Equation (|7]i, applied to the Chern-Simons action ( fTSl l, 
yields the (4 + 1)D analogue of Hall's law. 



u"-5pep(5)p(5) 

327r2 P' ■ 



(19) 



This equation, together with the conservation of the total cur- 
rent, = ibuik + .?'brane' rcproduccs the so-called chiral 
anomaly in (3 + l)D: a^j^mne = '^hE ■B,aH^ N/An^. 

Axion electrodynamics in (3 + 1)D can be recovered from 
the (4 + 1) -dimensional theory discussed here by dimensional 
reduction lfl3[[l4ll : Suppose that the five-dimensional electro- 
magnetic field is ^''-independent. Then /^(4) dx^A^, with 

7*^^^ a curve parallel to the z^-axis from one boundary compo- 
nent to the other one, plays the role of the axion field, in the 
(3 + 1) -dimensional action of axion QED, and the thickness, 
L, of the "slab" in (4 + 1) -dimensional space-time is related 
to the parameter I in front of the axion term in Eq. ( fTTb . The 
five-dimensional formulation shows that the time-derivative 
of the axion field plays the role of a (space-time dependent) 
chemical potential difference between left- and right-handed 
fermions |13, 14, 22]. 

Instabilities. The equations of motion derived from the 
Maxwell action and Eq. ( [TT] ), namely the Maxwell equations 
for the electromagnetic field and the equation 



la ^ ^ 



(20) 



exhibit an instability leading to the generation of helical mag- 
netic fields, as originally pointed out in Refs. lfl3l[l4ll : see also 
Ref. Ii22il . Linearization of the equations of motion around 
E = B = and a linear-in-time or time-periodic, but 
x-independent solution of Eq. (l20l l reveals unstable Fourier 
modes of the electromagnetic field, for small enough wave 
vectors (parametric resonance). Closely related instabilities 
have recently been discussed in Refs. [[32ii33[l . 

Conclusions. The key idea discussed in this letter is to 
promote global symmetries of systems of condensed mat- 
ter to local gauge symmetries and to study the response of 
such systems to turning on small, slowly varying gauge fields 
corresponding to those symmetries. By using general prin- 
ciples, in particular gauge invariance, anomaly cancellation, 
cluster properties and power counting, one is able to deter- 
mine the general form of the effective action or free energy 
of such systems in the scaling limit (as functionals of the 
gauge fields). This leads to a partial classification of states 
of condensed matter, including ones not characterizable by a 
local order parameter, in particular "topological phases". Of 
course, in applications to specific systems, the general con- 
siderations described in this letter must be supplemented by 
insights into structural properties that enable one to identify 
the relevant symmetries and associated gauge fields and the 
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physical meaning of the latter and to find out which terms in 
the general expression for the effective action (or free energy) 
are absent. 
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